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The aim of this paper is to investigate how the correlation properties of a stationary Markovian 
stochastic processes affect the First Passage Time distribution. First Passage Time issues are a 
classical topic in stochastic processes research. They also have relevant appHcations, for example, 
in many fields of finance such as the assessment of the default risk for firms' assets. 

By using some explicit examples, in this paper we will show that the tail of the First Passage Time 
distribution crucially depends on the correlation properties of the process and it is independent from 
its stationary distribution. When the process includes an infinite set of time-scales bounded from 
above, the FPTD shows tails modulated by some exponential decay. In the case when the process 
is power-law correlated the FPTD shows power-law tails 1/f'""'"^'''^ and therefore the moments (t") 
of the FPTD are finite only when n < (a — l)/2. We will also show that such power-law behaviour 
is not merely due to the fact that the process includes an infinite and unbounded set of time-scales. 
Rather, the time-scale must enter the FPTD with weights that must be distributed according to a 
power-law for large time-scales values. 

Finally, we will give a general result connecting the FPTD of an additive stochastic processes x{t) 
to the FPTD of a generic process y{t) related by a coordinate transformation y — f{x) to the first 
one. 

PACS numbers: 02.50.Ey, 05.10. Gg, 05.40.-a, 02.50.Ga 



I. INTRODUCTION 

Stochastic processes are used to model a great variety 
of systems in disciplines as disparate as physics [TH2 , ge- 
nomics [SI [3] , finance [101 IH] , climatology [12] and social 
sciences [T^ . 

Here we will devote our attention to the special class 
of Markovian stochastic processes that are also station- 
ary and that can be described by a Fokker-Planck (FP) 
equation [3]-[5]. Such processes are fully determined 
by the knowledge of their probability density function 
(pdf) W{x,t) and their conditional transition probabil- 
ity P{Xn+l,tn+i\Xn,tn). 

The simplest stationary Markovian stochastic process 
is the Ornstein-Uhlenbeck (OU) one that is characterized 
by an exponential autocorrelation function e^*/"^ where 
T is referred to as the time-scale of the process. For a gen- 
eral stationary Markovian additive process the method- 
ology of eigenfunction expansion [H [5] allows us to write 
the autocovariance function as the infinite weighted sum 
of exponential functions e~^* each characterized by a 
time-scale A^^, where the A values are the eigenvalues 
of the FP operator. By analogy with respect to the 
OU process, for an additive stationary Markovian pro- 
cess the time-scales are thus defined as the inverse of 
the eigenvalues of the FP operator. The methodology 
of eigenfunction expansion shows that the typology of 
time-scales included in a stochastic process is therefore 
strictly linked to the auto correlation properties of the 
process itself. The generalization to a multiplicative sta- 
tionary Markovian process has been considered in Ref. 
[I4l [TS] , where we considered stationary Markovian pro- 
cesses that include multiple time-scales [I6H2T] . In fact, 
we considered explicit examples of (?) short-range cor- 



related processes with an autocorrelation function that 
shows a decay modulated by some (stretched-) exponen- 
tial function, (ii) short-range correlated processes with 
an autocorrelation function that shows a power-law de- 
cay and (Hi) long-range correlated processes, thus show- 
ing anomalous diffusion 

By using such explicit examples of stationary Marko- 
vian stochastic processes, the investigations performed 
in Ref. [I4l [T5] mainly regarded the relationship between 
the stationary pdf or the conditional probability on one 
side and on the other side the typology of time-scales 
included in the process, i.e. its correlation properties. 
The aim of this paper is to investigate how the inclusion 
of more and more time-scales into stationary Markovian 
stochastic processes affects the First Passage Time dis- 
tribution (FPTD). 

First Passage Time issues are a classical topic in 
stochastic proceses research 27^ 2Sj • They also have rele- 
vant applications, for example, in many fields of finance. 
Just to mention, the assessment of the default risk for 
firms' assets is essentially based on the evaluation of the 
FPT under the assumption that the asset value V follows 
a geometric brownian motion |30H32j . i.e. an Ornstein- 
Uhlembeck process for the variable x = log(V^). There 
has therefore been much effort in obtaining the FPTD of 
the OU process thus neglecting the possible im- 

plications of the correlation properties of the asset value 
at the level of the FPTD and therefore at the level of the 
risk assessment. 

In this paper we will start considering the effects of 
the autocorrelation properties of a stationary Markovian 
stochastic process on the FPTD. Deliberately, our inves- 
tigations will be at a theoretical level, with the aim of 
understanding what are the relevant aspects of the prob- 
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lem. To this end we will concentrate on the First Passage 
Time through two barriers located at a; = a and x = b 
starting from position G [a, In fact, as explained 
below, such problem can be treated analytically. The 
search for a more realistic model in the context of de- 
fault risk or for other financial variables 29J will be left 
for the future. 

In this paper we will show that the large time depen- 
dance of the FPTD crucially depends on the correlation 
properties of the process and it is independent from its 
stationary distribution. Moreover, in the case when the 
process is power-law correlated the FPTD shows power- 
law tails As a result, in the case when the 
process is power-law correlated the moments (f") of the 
FPTD are finite only when n < (a — 1)/2. Finally, we will 
give a general result connecting the FPTD of an additive 
stochastic processes x{t) to the FPTD of a generic pro- 
cess y{t) related by a coordinate transformation y = f{x) 
to the first one. 

The paper is organized as follows. In section |Tl] we 
will give explicit examples of stationary Markovian pro- 
cesses that include multiple time-scales. We will consider 
short-range correlated processes that include an infinite 
set of time-scales bounded from above and processes that 
include an infinite set of time-scales that extends up to 
infinity. For both cases we will consider the case when the 
stationary pdf is gaussian or with power-law tails. In sec- 



tion III we will consider the mean FPT for the processes 
considered in section |ll] We will show how the simple 
knowledge of the mean FPT may be not enough to let 
us discern between processes with different correlation 



properties. In section IV we will derive a general expres- 
sion for the FPTD of stationary Markovian stochastic 
processes. By using the explicit examples of section |TT] 
we will then show how the correlation properties of the 
considered process affect the tails of the FPTD. Finally, 
in section |V] we will draw our conclusions. 



II. PROCESSES 

In this section we will briefly illustrate the stochastic 
processes to be considered in this paper. Further details 
can be found in Ref. [13] and Ref. [15] . 

Let us consider stationary Markovian processes that 
can be described by a FP equation [3]: 

^ W{x, t) = -^{h{x)W{x, t)) + ^{g{xrW{x, <)) (1) 



dt 



dx^ 



where h{x) and g{x) are the drift and diffusion coefficient 
appearing in the associated Langevin equation: 



x{t) = h{x) + 9{x)T{t) 



(2) 



where r(t) is a ^-correlated Gaussian noise term. 

For a stochastic process with a constant diffusion co- 
efficient, the eigenvalue spectrum of the FP equation de- 
scribing a stationary process consists of a discrete part 



Ao = 0, Ai,..., Ap and a continuous part ]Ac,-l-oo[ (Ac ^ 
Xp) associated with eigenfunctions The FP equation 
with constant diffusion coefficient can be transformed 
into a Schrodinger equation [¥ with a quantum potential 
Vs{x) = h{xy /A + dxh{x) /2. The eigenvalue spectrum 
of the Schrodinger equation is equal to the eigenvalue 
spectrum of the FP equation. The relation between the 
eigenfunctions of the FP equation and the eigenfunctions 
■0A of the Schrodinger equationis ipx = V'a'00- 

For a stationary process such that Vs{x) only admits 
a null eigenvalue and a continuum part of the spectrum, 
one can show that: 



{x{t + T)x{t)) - {x{t)){x{t + T)) = 



+ 00 



„2 -X 



^dX, (3) 



where cx = dx x (px (x) . This shows that the auto- 
covariance function is the infinite weighted sum of expo- 
nential functions e~^* each characterized by a time-scale 
A~^, where the A values are the eigenvalues of the FP 
operator. 

Eq. ([3]) clearly shows that the typology of time-scales 
included in a stochastic process is strictly linked to the 
auto correlation properties of the process itself. A gen- 
eralization to a multiplicative stationary Markovian pro- 
cess has been considered in Ref. [HI dS] • 



A. A short-range correlated process 

Let us consider the stochastic process described by the 
following Langevin equation [4]: 



i{t) = ~h{x{t)) +Dr{t), 
-fc if a; < 0, 

h{x) = ' 

-fc if a; > 0, 



(4) 



where is a real constant and T{t) is a (5-correlated gaus- 
sian noise term. The diffusion coefficient D will be set 
to unity hereafter. It is straightforward to show that the 
stationary distribution of the process is: 



Wsix) 



^exp(-^|a;|). 



(5) 



The above process is a stochastic Markov process de- 
scribed by a FP equation with constant diffusion coeffi- 
cient. As such, the FP equation maps onto an equiva- 
lent Schrodinger equation [4 with an appropriate quan- 
tum potential Vs{x). In this case the potential is given 
by: Vs{x)— — kS{x). It admits a ground state 

with null eigenvalue and a continuum set of eigenvalues 
A > /c^/4 whose eigenfunctions are given in Ref. [4]. By 
using the methodology of eigenfunction expansion [H [S] 
it is possible to prove that the autocovariance function 
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Rs{t) = {x{t)x{0)) of the above process is: 



= - 2^ + + -r3) (l - El•f(^/7))+ (6) 



A;2 



(2t- 1)(3 + 2t) exp(-r), 



which behaves hke a power-law with an exponential trun- 
cation for large time lags: Rs{t) ~ exp(— T)r^'^/2 as 
t — oo. As a result, this is a short-range correlated 
process characterized by an infinite set of time-scales 
bounded from above, i.e. A"^ < 4/^^. 



B. A power-law correlated process 

Let us consider the stochastic process described by the 
following nonlinear Langevin equation |141 120| : 

x = h{x{t)) + r{t), 

( -2^A^tan(v^x) if \x\ < L, 

h{x)^{ , (7) 

[ {1- V'^ + 4Vi)/x if |a;| > L, 

Vi L^/Votan(^/VoL) (l + LV^tan(VK)^ 



where L and Vq are real arbitrary constants and T{t) is 
a (S-correlated gaussian noise term. It is straightforward 
to show that the stationary distribution of this process 
is: 



Wiix) 



Nil cos{y/VQx) if \x\<L, 



(8) 



Nil/x°' if \x\ > L, 



where Nj and Njj are real constants that can be obtained 
by imposing that Wi{x) is continuous and normalized to 
unity. 

The quantum potential associated with this process is 
given by: 



Vsix) 



-Vo if < L, 
Vi/x^ if |a;| > i. 



(9) 



Such potential admits a ground state with null eigenvalue 
and a continuum set of eigenvalues A > attached to the 
null eigenvalue. 

The eigenfunction of the ground state is ipo = 
Bq cos(v^a;) whereas for \x\ > L it decays according to 
■00 = Aq x~"/^. The constants Aq and Bq are set by im- 
posing that TpQ is normalized and continuous in a; = ±L. 
It is worth noting that for > L the stationary pdf 
W{x) of the stochastic process is a power-law function 
decaying as |a:|~". The normalizability of the eigenfunc- 
tion of the ground state is ensured if a > 1. In the 
present study we consider stochastic processes with finite 



variance which implies a > 3. For |a;| > L the eigenfunc- 
tion ipx is a linear combination of Bessel functions ipx = 
Ax ^/xJ^{y/Xx) + Bx ^/xYt,{\/\x) where v = {a + l)/2. 
For |a;| ^ L wc find ■0^°'^'''' = Dx sin(-\/Vo + Xx) and 
^(^eiien) _ cos(\/Vo -l- X x) . The Coefficients Ax, Bx, 
Dx and -Fx are fixed by imposing that tpx and its first 
derivative are continuous in x = ±L and that ipx are 
orthonormalized with a ^-function of the energy. 

By using the methodology of eigenfunction expansion 
[U [5] it is possible to prove that the autocovariance func- 
tion {x{t + T)x{t)) oc T~^, where P = [a — 3)/2 thus 
showing that we are dealing with a power-law correlated 
stochastic process. In the range 3 < a < 5 the process 
is long-range correlated. As a result, this is a power-law 
correlated process admitting an infinite and unbounded 
set of time-scales. 



C. A short-range correlated process with a 
gaussian pdf 



Let us now consider the process of section II A and the 
coordinate transformation 1151: 



x^y^ fg{x) = %/2^Erf"^ [l - e^*^!^!] 



(10) 



By using the Ito lemma, one can show that, starting from 
the process of Eq. Q, in the coordinate space y = fg{x) 
one gets a multiplicative stochastic process whose sta- 
tionary pdf is given by: 



WM = 



1 ^ 1 2^ 

/27rs 2s 



(11) 



with s an additional arbitrary parameter. The corre- 
sponding drift and diffusion coefficients are given by: 

H{y)^^k'ey'/'^'^^ (l-Erf(H/V2;)) x (12) 

(l - Erf (|y|/\/2^) ) -e 2\/27rs) 



G(y) 



X I TT e 
ITT 



k ^/s e 



1-Erf |t/|/\/2; 



(13) 



where e = -|-1 when y > and e = —1 when y < 0. For 
large values of y the drift and diffusion coefhcient behave 
as: 



H{y) cx —k^s 



1 



G{y) 



-k. 



1 



y 



+00. 



According to the eigenfunction expansion methodol- 
ogy [H [5], the autocorrelation function of the process 
defined by Eq. Q with the coordinate transforma- 
tion of Eq. ^ is given by pg(r) = {{y{t + T)y(t)) - 
{yim/iiy'm- {yim where: 



{y{t)y{t + T)) = 



dXCle 



Ac 



Ac = fcV4, 



Cx = 



A": 

+00 



(14) 



dxfgix) iJo{x)^pxix), 
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where V'o(a^) and are the eigenfunctions of the 

Schrodinger equation with potential Vs{x)= k'^/A—k S{x) 
associated to the stochastic process of Eq. Q. 

As explained in Ref. [TS] such process admits 
an infinite set of time-scales bounded from above 
and its autocorrelation function decays like Rg{t) ~ 
exp(— fc^/4i)t~'^/^ for large time lags. As a result, this 
is a gaussian distributed, short-range correlated process, 
characterized by an infinite set of time-scales bounded 
from above. 



D. A short-range correlated process with a 
power-law decaying pdf 

Let us now consider again the process described by Eq. 
Q . By performing the following coordinate transforma- 
tion: 

(^^a)l/(a-l)gfex/(a-l)^ \x\ > L / D , 

x^y = fp{x) = \ (15) 
Dx, \x\ < L/D, 

and by using the Ito lemma, one can show that in the 
coordinate space y = fp{x) one gets a multiplicative 
stochastic process whose stationary pdf is given by: 



E. 



kL° -kL/D J, 
2D ^ y- 



e 7^, \y\ > L, 



(16) 



k ^-ky/D 



2D 



\y\<L, 



with L and a positive real constants. The requirement 
that the diffusion coefficient is continuous in y = ±L 
gives the following constraint: 



k = {a-l) 



D 



(17) 



Correspondingly, the drift and diffusion coefficients are 
given by: 

" -{a-2)/L^D^y, |y| > L, 
H{y) = { (18) 
+ {a-l)/LD^ |y|<L, 

D/L\yl |y| >L, 
G{y) - <; (19) 
D, \y\<L, 

According to the eigenfunction expansion methodol- 
ogy [H |5j, the autocorrelation function of the process 
defined by Eq. Q with the coordinate transforma- 



tion of Eq. (15 1 is given by Eq. (14) with fg{x) re- 



placed by fp{x). The relevant integrations can be per- 
formed analytically. One can therefore show that such 
process admits an infinite set of time-scales bounded 
from above and its autocorrelation function decays like 
Rp{t) « exp(— fc^/4i)t~'^/^ for large time lags. As a re- 
sult, this is a short-range correlated process with a power- 
law decaying stationary pdf, characterized by an infinite 
set of time-scales bounded from above. 



A power-law correlated process with a gaussian 
pdf 



Let us now consider the process of section [IlB| and the 
coordinate transformation (15j : 



x^y^ fi{x) 
r{x) 



p{x) 



a-1 ' 
2NiL 



// = \/2^Erf"'[l-r(x)], \x\ > L, 

fii^V2s^if"[l-p{x)], \x\<L, 

(20) 

{x-L) + 



Ni 



COs(yV^i)2 

Ni 



cos(^/T^i)2 



(sin(2yi^x) - sin{2^/VoL)y 



By using the Ito lemma, one can show that, starting from 
the process of Eq. ([t]), in the coordinate space y = fi{x) 
one gets a multiplicative stochastic process whose sta- 
tionary pdf is given by Eq. ( 11 ), with s as an additional 
arbitrary parameter. One can easily show that G(y) is 
continuos along the real axis, although its first derivative 
is discontinuos my — ±/(L). The drift coefficient H{y) 
is discontinuos in y — ±f{L). For large values of y the 
drift and diffusion coefficient behave as: 



1 



y 



-H"(y) « ^;+rexp( 



G(y)cx 



exp( 



'2s{a- I)' 



y +O0. 



According to the eigenfunction expansion methodology 
m \5\ , the autocorrelation function of the process defined 
by Eq. ^ with the coordinate transformation of Eq. 
(20) is given by Eq. (14) with fg{x) replaced by fi{x), 
Ac = and the eigenfunctions ipQ and ipx are solutions of 
the Schrodinger equation with the quantum potential of 
Eq. ^. 

As explained in Ref. [H] and Ref. [15] such process ad- 
mits an infinite and unbounded set of time-scales and its 
autocorrelation function decays like Ri{t) sa l/t^"^^)/^ 
with logarithmic corrections for large time lags. As a re- 
sult, this is a gaussian distributed, power-law correlated 
process, characterized by an infinite and unbounded set 
of time-scales. 



III. THE MEAN FIRST PASSAGE TIME 

In order to understand how the correlation properties 
of a stochastic process affect its FPTD we start investi- 
gating the mean FPT, which is the first moment of the 
FPTD. We will show in this section that processes with 
different correlation properties might show a similar be- 
haviour in the mean FPT. 

Let us consider the mean time T3.(A) that is needed 
to reach for the first time position a; ± A starting from 
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position X. This is obtained by solving the equation [S]: Eq. (271 and Eq. (24) clearly show that for large values 



h{x) 



dx 



dx^ 



(21) 



with boundary conditions Tx{A) = when x = ±A. Here 
h{x) and g{x) are the drift and diffusion coefficient ap- 
pearing in the non-linear Langevin equation associated 
to the considered stochastic process. 



A. Explicit examples of processes with a gaussian 
stationary distribution 

In this section we will compare the mean FPT of the 
two stochastic processes with gaussian pdf considered in 
the previous section. Indeed, the results shown here have 
already been considered in Ref. [Ej. However, for the 
sake of clarity we present them here again in a shortened 
form. 

For the stochastic gaussian process of section [TTC the 
mean FPT can be analytically computed by using the 
drift coefficient of Eq. ( 12 ) and the diffusion coefficient 
of Eq. (13) into Eq. (21). For the case when a; = one 
gets: 



To (A) 



1 

fe2 



log(^l -Erf(A/%/2 
1 



(22) 



1 



1 - Erf (A/A/2i) 



An alternative way to obtain such result is to analytically 
compute the mean FPT for the stochastic process of Eq. 
^ and then performing the substitution: 



1 



Ah-> --log(l-Erf(A/V2^) 



(23) 



For large values of A the mean FPT of Eq. ( 22 ) gives: 

A 



MA) 



z e 



/2s 



(24) 



For the process of section [II E| the mean FPT can not 
be computed by analytically solving Eq. (21 1. How- 



ever, one can follow the alternative way mentioned above. 
Specifically, we first computed analytically the mean first 



passage time of the process of Eq. (|7[) by using Eq. (21 1 



Subsequently, according to Eq. ( 20 ) , we performed the 
substitution: 



A^/7'(A), 
A^/7/(A), 



if A > L, 
if A < L. 



(25) 
(26) 



In the region A > fi{L) it is possible to obtain an ana- 
lytical expression showing that for large values of A one 
gets: 



T^{A)^[ze^y-\ 



II E has a mean FPT larger 



II C| This confirms that the 



of A the process of section 
than the process of section 
process that incorporates more time-scales is the slowest. 

In the two examples considered here, the way the cor- 
relation properties of the process affect the mean FPT is 
clear. In the case of Eq. (24) we have a linear growth in 

terms of the 



function, while in the case of Eq. \21\ 
we have a power-law growth with exponent larger than 
unity. 



B. Explicit examples of processes with a power-law 
decaying stationary distribution 

Let us now consider the mean FPT of the two stochas- 
tic processes with power-law decaying pdf considered in 
section HIl 

For the stochastic process of setion [ITD] the mean first 
passage time can be analytically computed by using the 
drift coefficient and the diffusion coefficient of Eq. ( 18 ) 
and Eq. (19 1 into Eq. (21 1. It is possible to obtain an 



analytical expression showing that for large values of A 
one gets: 



To (A) cx A"- 



(28) 



where a is defined in Eq. (17). For small values of A one 
gets To(A) « AV2. 

As mentioned in the previous section, for the stochastic 
process of section [TlB| the mean FPT can be analytically 
computed by using the drift coefficient and the diffusion 
coefficient of Eq. ([7| into Eq. (21). It is possible to 



obtain an analytical expression showing that for large 
values of A one gets: 



ro(A)aA 



Q + l 



(29) 



For small values of A one gets To(A) sa A^/2. 

Both in Eq. ( 28 ) and ( 29 ) the parameter a is the expo- 



nent of the pdf tail. Therefore, as for the processes with 
gaussian pdf, the mean FPT of the process that includes 
more time-scales is larger than the other. However, in 
the two cases considered here a power-law dependance 
from A is found. The qualitative behaviour of the mean 
FPT in both processes is therefore very similar, although 
the correlation properties of the two processes are very 
different. In fact, we will show that the different corre- 
lation properties of the stochastic processes display their 
effect at the level of the higher order FPTD moments. 



IV. THE FIRST PASSAGE TIME 
DISTRIBUTION 

A. A Theoretical result 

Let us suppose to be interested in the distribution 
'Da.b{xQ,t) of the first times at which a process reaches 
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the absorbing barriers x = a or x = b starting from a 
generic position xq with a < xq < b. According to Ref. 
[5] one gets: 

Ga,b{xo,t)^ f dxP{x,t\xo,0), (30) 

J a 

where P{x, t\xo, 0) is the conditional probability of find- 
ing the process in position x at time t, knowing that 
it was in position xq at time i = 0. According to the 
methodology of eigenfunction expansion 4, ,5 , for an ad- 
ditive stationary Markovian process described by a FPE 
admitting only a continuum set of eigenvalues, the con- 
ditional probability is given by: 



Schrodinger equation with the potential Vs{x) associated 
to the additive process x{t). As a result, the FPTD of 
the transformed process y{t) is given by: 



P{x,t\xo,0) = Mxf 



V'o(a^o) 



where the A-integration is extended over the continuum 
part of the spectrum and tpo{x), are the eigen- 

functions of the quantum potential Vs associated to the 
stochastic process considered. If the potential Vs{x) ad- 
mits a discrete part of the spectrum then an additional 
term involving a sum over all the discrete eigenvalues 
should be considered. By using the above expression for 
the conditional probability one finally gets: 



1 



-00(2^0) Jx 

X / dxil;Q{x)'4>x{x). 



dXe Tpxixo) X 



(32) 



The role played by each time-scale in the FPTD can be 
emphasized by rewriting Eq. ( 32 ) in the form 



Ga,bixo,t) 



dXe~^*yVx{xQ;a,b), 



yVx{xo;a,b) = ^-^^^"^ / dx il;o{x)tpx{x), (33) 
iPo{xo) J a 

where yVx{xo',a,b) can be considered as the weight with 
which each time-scale enters the FPTD. 

The above formulae can be used to obtain the FPTD 
in the case of the additive stochastic processes of section 
IHAland section HHl 

Let us now consider a stochastic process y(t) obtained 
by performing a coordinate transformation y = f{x) 
starting from a stochastic additive process x(<), as for 
the processes of section |II C[ section [TlD| and section [TlE| 
The conditional probability is given by 15 : 



P{y,t\yo,Q)^W{y) + 
1 i^o{x) 



(34) 



-At 



^ , ( N / d\^x{x)ipx{xQ)e , 

x = r^{y), XQ^f^'^{yo), 

where W{y) is the stationary distribution of the pro- 
cess y{t) and iPq{x), ipxix) are the eigenfunctions of the 



1 d 



dXe ipx{xo) X 



-00(3:0) dt Jx, 

X / dx'ipoix)'>pxix), (35) 



Xo = f ^ivo)- 



A direct comparison of Eq. ( 35 ) and Eq. ( 32 1 shows 
that: 

Kbiyo^*)='Df~Ha)j-Hb){r\yo),t). (36) 

Therefore, for a general stochastic process y{t), the prob- 
lem of finding the FPTD can always be re-conducted 
to the problem of finding the FPTD of an appropriate 
additive process. When the additive process admits a 
quantum potential Vs{x) that is exactly solvable then 
Eq. (33) can be used. If Vs{x) is not exactly solvable. 



one might either take advantage of the approximations 
schemes widely used in quantum mechanics to estimate 
the eigenfunctions or revert to numerical simulations of 
the process to estimate the FPTD. 



B. Explicit examples 

Let us now consider the stochastic process of section 
Moreover, let us (lx Xq = and b = A, a = —A 
in section |III| The integration over x required in 

It gives that 
Therefore 



II A 



as 
Eq, 



p2| can be performed analytically. 

fc2/4 when A P/A. 



>Va(0;-A,A) « ^ 
the large time behaviour of the FPTD is given by [36 



3/2 



t — ?► 00. 



(37) 



According to the result of Eq. ( |36[ ), the same asymptotic 
behaviour is shared by the FPTD of the process of section 



lie and section II D although they have different pdfs. 

Let us now consider the stochastic process of section 
II B| for the case when xq = and b ~ A, a ~ —A. The 
integration over x required in Eq. p2[ ) can be performed 
analytically. One therefore obtains that Wa(0; — A, A) « 
^(a-3)/2_ Therefore, the asymptotic behaviour of the 
FPTD is given by: 



V. 



A,A 



(0,0 - iA("+i)/2^ t 



(38) 



According to the result of Eq. (36), the same asymptotic 



behaviour is shared by the FPTD of the process of section 
|II E[ although they have different pdfs. In Fig. [TJwe show 



the function G_A.A(0,t) for the process of section II B 
The curves are obtained by using Eq. ([33]). We have 
considered the following three set of parameters: (i) L — 
1 and Vo = 0.985 (solid line) corresponding to a = 3.05; 
(ii) L = 1 and Vq = 1.16 (dashed line) corresponding to 
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law growing mean FPT. However, in the case of section 
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IIP I all the higher moments of the FPTD exist, while in 
the case of section II B only the moments (t") such that 
n < (a — l)/2 are defined. 

Finally, the result of Eq. ( 38 ) shows that the tail of 



the FPTD for power-law correlated stationary Marko- 
vian processes is different from what is expected for other 
non-Markovian processes with anomalous diffusion. For 
instance, in Ref. [35j it is shown that for a Fractional 
Brownian motion the tails of the FPTD decay like f^^"^, 
where 7 e [0, 2] is defined as {X{t)'^) cx H and {X{tf) is 
the mean square displacement of the process. 



Another process with an infinite and 
unbounded set of timescales 



FIG. 1: The figure sliows the function G_a,a(0, for the 
process of section |II B| The curves are obtained by using Eq. 
( 33 1 . We have considered the following three set of param- 
eters: (i) L = 1 and Vo = 0.985 (solid line) corresponding 
to a = 3.05; (ii) L — 1 and Vo ~ 1.16 (dashed line) corre- 
sponding to a = 4.00; (iii) L = 1 and Vo = 1.42272 (dashed 
dotted line) corresponding to a = 6.00. For all the three sets 
we considered A = 2L. In the bottom part of the legend 
we report the result of a nonlinear fit performed in the range 
t £ [20, 100]. These results are clearly in agreement with those 
ofEq. pel). 



a = 4.00; (iu) L = 1 and Vb = 1.42272 (dashed dotted 
line) corresponding to a = 6.00. For all the three sets we 
considered A = 2L. In the bottom part of the legend we 
report the result of a nonlinear fit performed in the range 
t £ [20, 100]. These results are clearly in agreement with 
those of Eq. (ISSl). 



The above results explicitly show that the tail of the 
FPTD is independent from the specific stationary distri- 
bution of the stochastic process considered. Rather, it is 
strictly related to the correlation properties of the pro- 



cess. In fact the processes of section II B and section HE 
have the same FPTD tail although their pdf is different. 
Similarly for the three processes of section |II A[ section 
EC] and section HlDl 

Another relevant comment regards the moments of the 
FPTD distribution. The existence of power-law tails in 
the FPTD of power-law correlated stationary Markovian 
processes implies that the moments of the distribution 
(t") are defined only for n < (a — 1) /2. In the case of the 
processes of section |II A[ section II C| and section |IID[ 



the exponential cut-off of Eq. ( 37 ) ensures that the mo- 
ments of the FPTD are always finite. That marks a rel- 
evant difference between a process admitting a bounded 
set of timescales and a process in which all timescales 
are included. As shown in section IIIIBl such difference 



is not caught by looking at the mean FPT only. In fact 
both the processes of section |II D| and section |II B[ both 
characterized by a power-law decaying pdf, have a power- 



The fact that the process of section |II B| includes an 
infinite and unbounded set of time-scales is a necessary 
condition for the existence of power-law tails in its FPTD. 
This is ultimately due to the fact that the eigenvalues 
spectrum of the Schrodinger potential of Eq. ^ is given 
by a continuum set of eigenvalues A > attached to 
the null eigenvalue associated to the ground state and 
therefore it includes an infinite and unbounded set of 
time-scales.. One might therefore wander whether or not 
such power-law decay is merely due to the fact that the 
process incorporates an unbounded set of timescales or 
something else must be considered. 

In general, any Schrodinger potential that asymptot- 
ically decays like l/x^^ would give a stochastic process 
with an infinite and unbounded set of time-scales. In 
the case fi = 1, we will show here that the correspond- 
ing FPTD tail does not necessarily display a power-law 
decay. 

Let us now consider the additive stochastic process as- 
sociated to the Schrodinger potential |TS]: 



Vs{x) 



-Vo a \x\ < L, 
V,/\x\ ff \x\ > L, 



(39) 



where L, Vq and Vi are real positive constants. This 
quantum potential is exactly solvable. In the region x > 
L the eigenf unctions are given by: 



M^) ^ Ao Vx Ki{2 y/Vi Vx), 



(40) 



iy/Xx 



AA-Fi,i(l-i^ -^,2,2iVXx) + 



C\U{l-t^ -^,2,2iVXx))+c.c 
2 V A ^ 



where Ki{-) is the K-Bessel function of order 1, [/(•) is 
the confluent hypergeometric function and -F'i,i(-) is the 
Kummer hypergeometric function [37] . 

The normalization constant Aq is fixed by imposing 
that the ground state is normalized to unity. The other 
normalization constants and C\ are fixed by imposing 
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that the odd and even eigenf unctions fulfill the normal- 
ization condition J dxipx{x)'il'\'ix) = S{X — X'). The drift 
coefficient h{x) aymptotically decays like h{x) w l/^/x. 

As mentioned above, the large time behaviour of 
^'-a.+aIO, t) is determined by the small energy behaviour 
of yV\{0; — A, A) [36]. One can explitely show that in the 
limit when A is small then: 

^ V A 

Ca C/i,i(1 - ^,2, 2i^/Xx) w ^/^ Kl{2^/Vlx), 

2 V A 

where /i(-)and Ki{-) are Bessel functions of order 1. 
As a result, the small energy limit of tpx is essentially 
determined by the small energy limit of ^^-^1,1(1 — 
-^,2,2i^/Xx). By making use of Eq. (13.6.8) and 
Eq. (14.1.7) of Ref. [37] one gets: 



Wa(0;-A,A) - Ao^a(O) 



x(41) 



dy^K,{2^,y) Fo 



Vi 

2VX' 



'Ay), 



where Fq is the regular Coulomb function with null index. 
By making use of Eq. (14.4.1) of Ref. [37] one gets: 

Wa(0;-A,A) = AoFfMAVA(O) X 

/A 00 
dy^K^[2^/Viv)Y,hkIk{2^/Viy), 

where /fc(-) is the Bessel functions of order k and the 
coefficients are defined in Eq. (14.4.6) of Ref. [57] . 
We have thus factorized the A dependencies from the 
y dependencies. As a result, the large time behaviour 
of the FPTD is given by the small energy behaviour of 

A\ipx{0) ~ e^N/A. By using the results of Ref. [36] one 
finally gets: 



A,A 



(0,0 



K = 3 



Stt \ 2/3 



The above result shows that the mere existence of an infi- 
nite and unbounded set of timescales is not enough to en- 
sure an asymptotic power-law behaviour in 2?_A^A(0,i). 
In passing, by using the same arguments it is also pos- 
sible to prove that the stochastic process associated to 



the quantum potential of Eq. ( 39 ) is not power-law cor- 



related, even though it incorporates an infinite and un- 
bounded set of timescales [21] . 



CONCLUSIONS 



The aim of this paper was to investigate how the FPTD 
is affected by the typology of time-scales included in a 
stochastic process, i.e. by the correlation properties of 
the process itself. The motivations are both theoretical 
and in view of possible applications. In fact, First Pas- 
sage Time issues have relevant applications in many fields 
of finance such as the assessment of the default risk for 
firms' assets. 

In this paper we have shown that the tail of the FPTD 
crucially depends on the correlation properties of the pro- 
cess and it is independent from its stationary distribu- 
tion. When the process includes an infinite set of time- 
scales bounded from above the FPTD shows tails mod- 
ulated by some exponential decay. In the case when the 
process is power-law correlated {x{t)x{t+T)) '-^ 7-"'-""'^^/^ 
the FPTD shows power-law tails l/t^^+^^Z^. As a re- 
sult, in the case when the process is power-law corre- 
lated the moments {t") of the FPTD are finite only when 
n < (a — l)/2. We have also shown that such power-law 
behaviour is not merely due to the fact that the process 
includes an infinite and unbounded set of time-scales. 
Rather, it is also required that each time-scale enters the 
FPTD in a peculiar way. In fact, in order to observe 
a power-law tail in the FPTD the weights 'yV\{xQ;a,b) 
must be distributed according to a power-law for small 
values of A, i.e. for large time-scales. 

Moreover, we have given a general result showing that 
given (i) and additive stochastic processes x{t) and (ii) 
another process y{t) related by a coordinate transforma- 
tion y = f{x) to the first one, the FPTD of y{t) can be 
obtained from the one of the x{t) process by using Eq. 
(36 1. Therefore, for a general stochastic process y{t), the 



problem of finding the FPTD can always be re-conducted 
to the problem of finding the FPTD of an appropriate 
additive process. When the additive process admits a 
quantum potential Vs{x) that is exactly solvable then 
Eq. (33) can be used. If Vs{x) is not exactly solvable. 



one might either take advantage of the approximations 
schemes widely used in quantum mechanics to estimate 
the eigenfunctions or revert to numerical simulations of 
the process to estimate the FPTD. 



Finally, the result of Eq. ( 38 ) shows that the tail of 



the FPTD for power-law correlated stationary Marko- 
vian processes is different from what is expected for other 
non-Markovian processes with anomalous diffusion, such 
as the Fractional Brownian motion. This result might be 
used to test the Markovian property of a generic stochas- 
tic process. 
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